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Semileptonic B decays into even parity charmed mesons 

Manuela De Vito and Pictro SantorcUi 
Dipartimento di Scienze Fisiche, Universita di Napoli "Federico H", Italy, 
Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Italy 

By using a constituent quark model we compute the form factors relevant to semileptonic transitions 
of B mesons into low-lying p-wave charmed mesons. We evaluate the dependence of these form 
factors and compare them with other model calculations. The Isgur-Wise functions ri/2 and 
are also obtained in the heavy quark limit of our results. 
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I. INTRODUCTION 

Recently, BaBar Collaboration has discovered a narrow state with = 0+ with a mass of 2317 MeV, .D*o(2317) 
0. The existence of a second narrow resonance, Dsj(2460) with = 1+, was confirmed by CLEO Ql- Both 
states have been confirmed by BELLE 3]. Soon after the discovery, another set of charmed mesons, Dq'^(2308) and 
Df (24:27) which have the same quantum numbers = (0"'", l"*") as Dgj has been discovered by BELLE j^. Before 
their discovery, quark model and lattice calculations predicted that the masses of these states, in particular D*o(2317) 
and D'gi{24:60) , would be significantly higher than observed |^,[^. Moreover, these states were predicted to be broad 
due to the fact that they can decay into D K and D* K, respectively. Experimentally, the masses of 15*0(2317) and 
D'gi{2A60) are below the DK and D* K thresholds and hence they are very narrow. These facts inspired a lot of 
theorists to explain the puzzle 

In this paper we will focus our attention on the weak semileptonic transitions of B mesons into lower lying n-wave 
charmed mesons (D**). These transitions were studied, within a quark model approach, for the first time in Q and, 
more recently, in Q where the authors take into account the symmetries of QCD for heavy quarks already used 
in The light-front covariant model was adopted to study the same subject in ^^l- The relevant form factors 
were also evaluated, in the framework of QCD Sum Rules J3, in 

Here we employ a simple constituent quark model Il7| to evaluate semileptonic form factors of B mesons into 
p-wave charmed mesons. The plan of the paper is the following. In the next section we describe our quark model; 
the third section is devoted to introduce and evaluate the s-wave to p-wave form factors. Our way to fix the free 
parameters of the model and the resulting form factors are discussed in section four, while in section five the heavy 
quark limit of the form factors are computed and compared with Heavy Quark Effective Theory predictions; the T1/2, 
and T3/2 are also evaluated. In the last section we show and discuss our numerical results. 

II. A CONSTITUENT QUARK MODEL 
In our model |l6lll7| any heavy meson H(Qq), with Q G {6, c} and q £ {u, d, s}, is described by the matrix 



V3' 



2mc 



2mo 



(1) 



where mg (m^) stands for the heavy (light) quark mass; , q2 are their 4— momenta (cfr Fig.QJ. 'ipnik) indicates the 
meson's wave function which is fixed by using a phenomenological approach. The meson constituent quarks vertexes, 
F in Eq. QJ, are fixed by using the correct transformation properties under C, P and to enforce the relation 



^|V'ff(fc)P-2Mff. 
For the odd parity, s-wave heavy mesons — (0^, 1^), F is given by 

To- = -«75i 
ri-(e, 91,92) = 



(2) 



2mQTOg 



LI LI 

gl - g2 

ruH + mq + ruq 



2 mQiriq 



(3) 
(4) 
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FIG. 1: Quark model diagram for the hadronic amplitude relevant to semileptonic M M' decay involving Qi — » Qa 
transitions. The thin lines represent quarks, the thick ones mesons. The gray disks represent the quark-quark-meson vertexes. 



where e is the polarization 4— vector of the (vector) meson H. The vertexes of the lower lying even parity heavy 
mesons, instead, are given by the matrices ^ 



0+ 



T3p^{e,qi,q2) 
ripi(£, 91,92) 



(mg - TO^)(gf - g^) 



3mnm„ 



75 A 



(gf - g2 ) "^g 



751 



6 mqinrLq 



(5) 
(6) 
(7) 



As already discussed in |l6l Il7j| . the 4— momentum conservation in the meson-constituent quarks vertexes can be 
obtained defining a heavy running quark mass. For details we address the reader to the references Here, for 

the sake of utility, we recall all the remaining rules of our model for the evaluation of the hadronic matrix elements 
of weak currents: 



a) for each quark loop with 4-momcntum k we have 



(fik 
(2^ 



(8) 



a colour factor of 3 and a trace over Dirac matrices; 
b) for the weak hadronic current, ^2 qi, one puts the factor 



/mi /ma 
El V E2 



(9) 



where with we indicate a combination of Dirac matrices. 



III. FORM FACTORS 



In this section we evaluate the form factors parameterizing the ^(0,1) and (0+, '^Pi, ^Pi) weak 
transitions. The decomposition of these matrix elements of weak currents in terms of form factors are the following 



^ We don't consider in this paper the tensor mesons. 
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(see also 



<h;,^{p')\v^\Ho-{p)> = 

<H[-ip\£)\V^-A^\Ho-{p) > = 

<HoAp')\A^\Ho-{p) > - 

<Hsp^ip',e)\V^~A^\Ho-ip) > - 

<Hip^ip',£)\V^~A^\Ho-ip) > = 



f+{q')ip^+p',.) + f-{q'){Pf.~p'^) , (10) 

2g(q2)eM-/3 elp^p'^ 

+i£*-p) [a+iq^)ip,+p'^)+a^iq^){p^-p'^)\} , (11) 
F+{q'){p,,+p'^,) + F^{q^){p,-p\,) , (12) 
~^{F'iq^)e; +is*-p) [A'^{q'){p,+p'^) + A'_{q'){p,-p',)]} 
+2G'iq') ef^^^f" slp^p'p, (13) 
-^{F{q')e; +{e*-p) [A+{q'){p,+p'^)+A.{q'){p,-p'^)]} 
+2 G{q^) e^--"^ . (14) 



^uPaPfi ■ 

The calculation of the form factors in Eqs. (|l()|l and I|1HI for the case of i? ^ D{D*) transitions has been done in 
Ref. 01 . However, for the sake of utility, the analytical expressions are reported in appendix 1X1 
One of the main results of this paper is the calculation of the form factors appearing in Eqs. 112(1 . (|13l) and (|14|l . By 
way of an example, in the following we describe the calculation of the matrix element < Hq+ {p')\Afj,\HQ- (p) > and 

give the expressions of the form factors F±. In appendix IbI we collect the expressions for G('\ F(') and A^'. Note 

that all the calculations are done in the frame where q^ = {E^, k ~ q), = {Ei, k), q!^ — {E2, ~k). 
Let us start considering the 0~ — > 0^ transition, 



<H„+{p')\Q3^^^5Qi\Ho-{p)> = - 

Tr 



-^V(^)V'o-(fc)W^^ 



(2^)3 
- ^2 



TO2 



(7oro+(g3,g2)'7o)^7Z — Tm75 



2m2 



2m3 



2m 



0-) 



2m2 



where V is the integration domain (see Refs. 0,^3) defined by 

Max(0,A:_)< k <Mm{KM,k+) 

Max(-l,/(fc,|g|)) < cos{e) < +1 

< < 27r 

with 



fc± - 
fikAq\) = 



|g| {nip + m^) ± {rrip - ml)^ 



2m], 



2y/m 



(m|, + m\) 



2k\q\ 



(15) 

(16) 

(17) 
(18) 



(j) and 6 are the azimuthal and the polar angles respectively for the tri-momentum k. Km = [fn] — ra^) / (2rai) and 
TO/ {mp) is the mass of the initial (final) meson: in Eq. (|15|) mp = TOq+. We choose the z— axis along the direction 
of q, the (tri-)momentum of the W boson (cfr Fig. P). 

The analytical expressions for the form factors can be obtained by comparing Eq. (|15|l with Eq. p2|) 



F+{q^) 



Pdkdcos9^Ji{k)iP*p{k) 



8 TO/ ^ {-dl2^ + TO/2) (tof^ - 323^) El 



E-i 



di2S2zE2 - mi (((ii3 - 2m2) m2 + to/£'2) + -—r cos{9) (-mp'^mj + di2S23 (to/ - Ep) + mj'^Ep) 



k1 



(19) 



F-il' 



where {di 



k'^dkdcosBiJji {k)ipp{k) 



8 7r2 TO/ ^ {-di2^ + TO/2) {mp^ - 523^) El E3 

k 

™2TO/Si3 + {mi^ + di2S23) E2 - cos(0) {mp'^mi + c?i2S23 (to/ + Ep) + mi'^Ep) 



(20) 



TOi 



TOj and 



TOi + TO j ) . 



The expressions for the remaining form factors in Eqs. pU|) - (|14|l are collected in appendix IXI and IbI 
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FIG. 2: The B — > D'Eft spectrum. The dashed (dotted) hne corresponds to the exponential (gaussian) vertex function. The 
data are taken from Ref. I2III 



TABLE I: The two best fit sets of values of the free parameters are obtained for the exponential (exp.) and gaussian (gauss.) 



Parameter 


fitted values (exp.) 


fitted values (gauss.) 


rriq 


23 MeV 


34 MeV 


LOB 


101 MeV 


108 MeV 


LUd = LUD' = LUD" 


51 MeV 


186 MeV 




0.041 


0.043 



IV. FIXING THE FREE PARAMETERS 



The numerical evaluation of the form factors given in Section lllll requires to specify the expression for the vertex 
functions and the values of the free parameters of the model. For the vertex functions we adopt two possible forms, 
the gaussian- type, extensively used in literature (see for example [T^ 'l 




^H{k)^A.jMH\l—^e^v\-7^\. (21) 



and the exponential one 



V'^(fc)=4^J^ exp(-A\ (22) 

which is able to fit the results of a relativistic quark model regarding the shape of the meson wave- functions [19| . 
In our approach ljh is a free parameter which should be fixed by comparing a set of experimental data with the 
predictions of the model. In this paper we choose to fix the free parameters by a fit to the experimental data on the 
Br{B — > Dii') j20j and on the spectrum of i? ^ D*£i' process .2 Ij . 

The quality of the agreement between fitted spectrum and the corresponding experimental data may be assessed 
by looking at the Figure [3 It should be also observed the very small differences between the B D*lv spectrum 
using the vertex functions in Eqs. (|21l) - (|22|l . Regarding the B Div branching ratio, we obtain 2.00(2.01)% for 
the exponential (gaussian) vertex function to be compared to the experimental value: Br{B Div) = 2.15 ± 
0.22 (2.12 ± 0.20)% for the charged (neutral) B meson. At this stage the two different form of the vertex functions 
agree equally well with the experimental data. However, differences emerge when single form factors are considered 
(cfr for example Table llTT|l . 
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V. HEAVY QUARK LIMIT 



In this section we perform the heavy quark hmit for the form factors obtained in the previous sections. Before to 
do this we briefly remind the imphcations of the HQET on the heavy meson spectrum. In the quark model, mesons 
are conventionally classified according to the eigenvalues of the observables J, L and S: any state is labelled with 
the symbol So, if we consider the lower lying even parity mesons (i = 1), the scalar and the tensor mesons 

correspond to '^Pq and states, respectively. Moreover, there are two states with J = 1: the ^Pi and ^Pi, they can 
mix each other if the constituent quark masses are different as in the case of charmed mesons. For heavy mesons the 
decoupling of the spin of the constituent heavy quark, sg, suggests to use a different set of observables: the total 
angular momentum of the light constituent, jq{— Sq + L), the orbital momentum of the light degree of freedoms 
respect to the heavy quark, L, the total angular momentum J (= jq + sq), any state is labelled with L'J . In this 

1/2 3/2 

representation the scalar and the tensor mesons are labelled with Pq and , respectively. The axial mesons are 
classified as Pj^^^ and P^"^ and are related to the "^Pi and ^Pi states by the 



^1 



3/2 



Pi> 



1/2 \ _ 



(23) 
(24) 



The scaling laws of the HQET concern the 0- ^ (Po'",P;'",Pi'") form factors, they are defined combining the 
form factors in Eqs. (O, (d) and the relations in Eqs. ^ and For example G'^''^{q^) = .JTJiG{q^) - 

To extract the heavy quark mass dependence from the expressions of the form factors wc follow the same approach 
used in our previous paper (17i |. We introduce the variable x, defined hy x — (2ak)/mF, in such a way, neglecting the 
light quark mass respect to the heavy ones, the integration domain, near the zero-recoil point, simplify to < x < a, 
< 9 < TT and < ^ < 27r. Therefore, if we look at the expressions of F±{q'^), Eqs. (|19|I - H2Q(I . near the zero recoil 



point (i. e. q"^ ~ Qmax) have, neglecting terms of the order of x^. 



F±{q') 

For a ^ 1 the integration can be easily done giving 

3/2 



dx 'ipi{k{x)) ippikix)) 



■mjp {iriF T ™/) (7 — 3 ui) x^ 
384 m/ TT^ 



(25) 



F±iq')\q^-^ql 



2\/2 ( 



1 — - (w — 1) ) + o((w — 1)^) gaussian — type 



exponential — type . 



(26) 



\{l-\{w~l)]+o{{w-lf) 



Thus, the Isgur-Wise function resulting from our model is given by 

11 19 

ri/2{w) = ^-^{^-l) + -^{w-l)'+ o{{w - If) , 

where we have also written the term {w — 1)^ which was neg lected in Eq. 



(27) 



A similar analysis can be performed on the heavy to heavy 
combination 



P^^^ form factors. Let us start considering the 



GV2(q2) 



1 



G{q') 



F'^Hq') F{q^) 



G'{q') 
F'{q') 



(28) 



which defines the ^ Pi form factors. It is very simple to obtain their scaling laws in the limit of heavy quark 
masses. Following the above method we obtain 



PV2(^2) 



N < 2^miraF{w — 1) 



(w- 1) +o{{w - 1)2), 



(29) 
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TABLE 11: The Isgur-Wise functions ri/2 and r3/2 at zero recoil and their slope parameters. 



^1/2(1) 


2 

P?/2 


^■3/2(1) 


2 

Pa/2 


Ref. 


0.33 


0.75 


0.83 


1.29 


This work 


0.34 


0.76 


0.59 


1.09 




0.22 


0.83 


0.54 


1.5 


[22] 


0.31 


1.18 


0.61 


1.73 


[23] 


0.13 ±0.04 


0.50 ± 0.05 


0.43 ± 0.09 


0.90 ± 0.05 


QCDSR [27] 


0.35 ±0.08 


2.5 ± 1.0 






QCDSR(NLO)f26] 


0.38 ±0.04 




0.53 ± 0.08 




Lattice [241 



where 



N 



2^2 f 



3/2 

gaussian — type 

(30) 

exponential — type . 



Similarly, we can evaluate the r3/2 Isgur-Wise function obtaining 

T3/2H = ^ - - 1) + - 1)2 + o{{w - \f) . (31) 

A comparison between our results and some others coming from quark models, QCD sum rules and Lattice calculations 
can be done looking at the Table ITll The values of the r functions at zero recoil point and their slopes are compatible. 
In particular, it should be observed that our results for ri/2 are practically the same obtained in the Isgur Scora 
Grinstein Wise (ISGW) model [1| and QCD Sum Rules findings [2^|23|.^ Regarding r3/2, our result at zero recoil 
point is slightly larger of the results coming from other models, while the slope is comparable with others. 

Relations between the slope of Isgur-Wise function and r functions at zero recoil points were derived, in the form 
of sum rules, by Bjorken [23 and Uraltsev 

„2 1 



P 4 



E 1-1/2(1)1' +2Ei-3;2(i)r (32) 

n 11 

\ - El-52(l)P-El-V2(l)^ (33) 

n n 

where n stands for the radial excitations and is the slope of the Isgur-Wise function i,{w) which, in our model p^ . 



m = 1 - ^(u; - 1) ± g(u; - \f ± o{[w \f) . (34) 

Our results for rt = oversaturate both the sum rules. For the Bjorken sum rule this is due to the small value we 
obtain for the derivative of the Isgur-Wise function (p^) which is in any case compatible with the experimental value 
— 0.95 ± 0.09 [23]. We plan to study this problem in a separate work. However, a detailed discussion on these sum 
rules and the findings of quark models can be found in [3nl |. 



VI. NUMERICAL RESULTS AND DISCUSSION 



All the results discussed in the previous section has been obtained without fixing the free parameters of the model. 
In this one we use the fitted values of the free parameters in Table (cfr section IIVI for discussion) to obtain the 



^ The values and the slopes of r functions are obtained fitting the numerical results obtained in Morenas et al. (see Ref Q) using the 
ISGW model. 
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TABLE III: B D** form factors evaluated at 5^ = and at q^ax = ("^s — m.D**)'^ by using the vertex function in Eq. 1221 . 
In parentheses the values obtained using the gaussian vertex function (cfr. Eq. H21|l ). 





This work 


Ref. \8] 


Ref. Ll^ 




f (U) 


\Hmax J 




^ \Hmax J 




P('n2 -) 

y^max J 


Fi 


-0.32 (-0.30) 


-0.35 (-0.33) 


-0.18 


-0.24 


-0.24 


-0.34 


Fo 


-0.32 (-0.30) 


-0.025 (-0.036) 


-0.18 


0.008 


-0.24 


-0.20 


/l(l/2) 


-0.25 (-0.23) 


-0.31 (-0.28) 


-0.18 


-0.39 


-0.075 


-0.083 


4(1/2) 


0.096 (0.088) 


-0.0018 (-0.0029) 


0.070 


-0.002 


0.073 


0.071 




0.69 (0.63) 


0.87 (0.79) 


0.49 


0.91 


0.32 


0.56 


y{i/2) 


0.67 (0.61) 


0.84 (0.76) 


0.44 


0.81 


0.31 


0.55 


4(3/2) 


-0.61 (-0.58) 


-0.81 (-0.77) 


-0.20 


-0.46 


-0.47 


-0.76 


^(3/2) 


-0.13 (-0.13) 


-0.016 (-0.023) 


-0.005 


-0.008 


-0.20 


-0.26 


^(3/2) 


0.70 (0.65) 


1.27 (1.19) 


0.33 


0.72 


0.25 


0.47 


y(3/2) 


-0.81 (-0.77) 


-1.09 (-1.03) 


-0.44 


-0.71 


-0.61 


-1.24 



TABLE IV: Parameters of the B — > D** form factors. The functional dependence is either polar : F{q^) — F(0)/(1 — 
aq'^ /m%) or linear: F{(f) — F(0)(1 -I- hq^). In parentheses the values obtained using the gaussian vertex function (cfr. Eq. 



Form Factor 


F(0) 


a 


b (GeV~2) 


Fi 


-0.32 (0.30) 


0.263 (0.233) 




Fa 


-0.32 (0.30) 




-0.109 (-0.103) 


4(1/2) 

^0 


-0.25 (-0.23) 


0.661 (0.626) 




^(1/2) 


0.10 (0.092) 




-0.120 (-0.120) 


yl(V2) 


0.69 (0.64) 


0.695 (0.680) 




^(1/2) 


0.67 (0.61) 


0.695 (0.679) 




4(3/2) 


-0.61 (-0.59) 


0.846 (0.834) 




^(3/2) 


-0.13 (-0.13) 




-0.102 (-0.0956) 




0.72 (0.67) 


1.53 (1.54) 




y(3/2) 


-0.81 (-0.77) 


0.872 (0.873) 





numerical results of the form factors. First of all, in Table Hill are collected the values of the form factors for the 
B D** transitions evaluated at zero recoil point {q^ax — {^b ^ nfiD—Y) and at = 0. We consider the charmed 
final state with the following masses: to(D5) = 2.40 GeV , m{D[) = 2.43 GeV, m{D{) = 2.42 GeV Note that 

we are considering, for a better comparison with other calculations, the elicity form factors (cfr for definitions, for 
example, ^Jl)- Looking at the Table ITTll we can see that the absolute values of our form factors (at q^ = 0) are larger 
than the ones in Ref |8lll3|. this naturally implies larger branching ratios in our model. In particular, our predictions 
on the branching ratios, using the exponential (gaussian) vertex function, are {tbo = 1.536 x 10^^^ s ^20] ) 



Br{B^ - 






= 2.3(2.1) X 10^3 


(|Kb|/0.041)2 


Br{B° - 


^ D'+r 




= 2.0(1.6) X 10-3( 


1^61/0.041)2 


Br{B° - 


D+t 




= 8.0(7.3) X 10-3( 


1^61/0.041)2 



Regarding the q^ dependance of the form factors, we find a very good agreement with numerical results assuming 



^ and D\ represent, respectively, the two different physical axial- vector charmed meson states. The physical D'^ ( Di) is primarily 

P\^^ (Pi^^)- They differ by a small amount from the mass eigenstates in the heavy quark limit, for a discussion see |3 J . In this paper 
we neglect these differences. 
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the following polar expression 



F(0) 



(36) 




the fitted values of a can be found in Table Hvl It is interesting to observe that the effective pole mass is not far from 
the mass of the Be meson. A different dependence exhibit the form factors Fq and Ai ; for them we use the form 



and the values of b are collected in Table HVI 

In conclusion we have obtained in a very simple constituent quark model all the semileptonic form factors relevant 
to the transition of B into the low-lying odd and even parity charmed mesons. The free parameters of the model 
have been fixed by comparing model predictions with the B D*£h' spectrum and B — > Dtv branching ratio. Our 
numerical results are generally larger than the results of other models. However, form factors reproduce the scaling 
laws dictated by the HQET in the limit of infinitely heavy quark masses. 




(37) 



APPENDIX A: FORM FACTORS /±, g, f, a± 



In this appendix we collect the analytical expressions for the 
Eqs. (fTn|l and lfTT|l . respectively. 



0" and 



1 form factors defined in 



Ml') 




X 



(Al) 



I 



dkdcosO 



X 




'23' - mp^) (^12' - mj2) El E3 {mp^ - Ep^) 




(A2) 
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^ sin{9y di2 E2 _ fc cos(e) (rf23 mi+dx2 Ep) 

!{q') ^ I edkdcosexPiikWpik)— '^^"^ , (A3) 



8 7r2 (mF + S23) y (^23^ - mj^^) (c^i2^ - 
{ — (to2'^ 7713) — fc^ TOi TOi? + mi mi?^ + TO3^ m/^ + 7713 m^? m/^ + nii nip^ S23— 

^23^ 'm-i {mp + S23) + ^23 (to_f + S23) + fc^ ™/ - m2^ (to_f ("t-3 + + -2 m/ i?_F)+ 

fc^ cos(2 0) (mi mp — m2 mp — mi Ep) + m2 (fc^ m/? + (ms — m/?) {m^ + mp)^— 

mp mi'^ + 2 {rri:} + mp) mi Ep)} , (A4) 



a+('7)+a-(g) = k^dkdcosO 



U2...„„„. ^7(fc)V'>(fc) 



47r2|9l^m/2 (m/^ + s23) (^23^ - mp'^) [di2^ - mj'^) Ei E3 
|-(fc^ mp^ {di2 mp - m/ Ep)) + \qf E2 (2 di2 mp E2 + m/ (-ms^ + (m2 + mp)^- 
2 i?2 Ep)) + k cos{9) (fc cos{d) {di2 mp — mj Ep) {mp'^ + 2Ep'^)+ 

\q\ Ep (—4 di2 mp E2 + mj {m^ ^~im2+mp)^+iE2Ep)))} , (A5) 
a+(g^) — a_ (q^) = / k^dkdcos6 



Mk)rp{k) 



Aii'^\q\^mi [mp + S23) \J (^23^ - mp"^) {di2^ - m/2) Ei E3 Ep 
|m/ (-(fc^ (|(7|^ + m_F^)) - (m2 {mp + S23) - mi E2)) Ep + di2 mp Qqf + mp'^)- 

\qf {mp + Sis) E2 Ep) + k cos{0) (fc cos(6l) (-(di2 mp) + mi Ep){\qf + mp'^ + 2Ep'^) + 

\q\ Ep {di2 (2 mp E2 + {mp + S13) Ep) + m/ (^23 {mp + S23) -{mi + 2E2) Ep)))} , (A6) 

where dij ^ mi — mj and = m^ + mj (with m^ the mass of z— quark); m/ and mp are the masses of the mitial 
and final mesons, respectively. Ep{= q"^ + m^) represents the energy of the final meson. The angle 9 is defined in 
Section Uni after Eq. ||THJ). 

APPENDIX B: FORM FACTORS G*'', F' ', A^-^^ 

In this appendix we give the expressions of the form factors appearing in Eq. H14() (0^ -^^Pi transitions). We use 
the same notations of the previous appendix. 

G{q') = - I k^dkdcose Mk)rp{k)Vik-sr.{e)'mp ^g^^ 

J 8^2 (_rf23' + mp^) J {-d,2^ + m,2) {mp^ - 823^) E^E^ 



F{e) . - / '-^ v--^ -f-, ( --) -/■: ( --) + ' + + "^^'^ (B2) 

A+(g2) + A_(,2) = - Ik'dkdcose Mk)rp{k)V^mp 



^E: 



8^2 l^f (_d232 + ^^2) mi^^{di2^ ~ mi^) {~mp^ + 523^) E^ 
|fc2sin(6')^mi^2 {-q^ + 2^12^23 + m^^ +m/2) - 2 ([(/l ^2 - A: C08{e) Ep) x 
(M ((^23^ - m^^) m/ + 2 di2 d23 E2 + {-q^ + mp"^ + m/^) E2) - 

k coa{e) {-q^ + 2 di2 d23 + mp^ + m/^) Ep) } (B3) 



A+{q^)-A.{q^) 
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fc2 dk dcos{0) {k)ip*p (fc) mp 



\^2miEF (\q\^ {di2^ - mi^) E2 + Ep (di2 ^23 + "^/ Ep) + 
k cos{6) (-3 k coa{9) Ep (di2 ^23 + mj Ep) + 

|gl {mi^Ep + 2miE2Ep + di2 (2 d23 E2 - Ep)))) | (B4) 
Regarding the form factors in Eq. (|13|l we have 



2n /3 f k dkdcos9'i/ji{k)ipp{k) 



{ - (M (('^23 (^^ + cJ23 "^2) - ^2 mp"^) mi + di2 (^23^ - mp^) E2)) + 

k coa{e) {k \q\ coa{e) ^23 w/ + (^23^ - mp'^) (mi S23 + ^12 -Ef)) | (B5) 



/3 /■ k^dkdcoseijiikUtik) 

F {q-} = 



8 ^2 (-^23' + 2) (^12' - mj2) (-^^.2 + ^232) E^E^ 
|-2 fc^ sin(0)^ ((ii2 mp"^ + ^23 "i/ -Bj") + (c^23^ - rap"^) (- (mi^ S23) + "^/^ S23 + ra\ [-mp'^ + 523^) + 
1112 {mp"^ - 11x3823 -'2 mi Ep))] (B6) 



3 [ k'^dkdcos0'ipi{k)'ijj*p{k) 



167r2|9f (-rf23'+mf^2)^^2,^(rf^22_^^2) {-mp^ + .s^^2) E,E3 

£'2 (m^^ TO/ S23 + 2 (ii2 mp'^ E2 + ^23 to/ (- ((^23 S23) + 2 £'2 -Ef)) - 
k \q\ cos{0) Ep [mp"^ mi S23 + 4di2 mp'^ E2 + ^23 to/ (- (^23 S23) +'^E2Ep)) + 
fc2 (di2TOF^ +d23TO/i;F) (-TOF^ +cos(6l)^ (tof^ + 2Sf^))) } (B7) 
A\{q^)~A' = + f edkdcos{e)iPi{k)rp{k) 



8 7r2 |g| (-^23^ + TOF^) TO/ E^F \l {di2^ - mi^) {-mp^ + 823^) E1E3 
|-(fc2 (di2 TOF^ + d23 TO/ £f) {-mp"^ + 2Ep'^)) + \q\^ (2 ^23^ TO2 to/ ^f + ^23 {k^mi+ 

2 (toi^ + TO2 TO3 — TO/2 — TOl S23) E2) Ep + TOf^ (fc^ dl2 — 2 (to2 To/ + dl2 E2) Ep)) + 

kEp {-3 k cos(2 6) Ep (o?i2 tof^ + ^23 TO/ Ep) — 

2 |g| cos(6') (c?23^ TO/ S23 - tof^ (to/ S23 + 2 di2 E2) + 

d23 mi' Ep + (to22 ms — toi (tof^ + 0^23 S23) + TO3 to/ (to/ + 2 E2) - 

m2 {ms^ -mp^ + mi^ + 2miE2)) Ep)) ^ (B8) 
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